Re normalization of pinned elastic systems: how does it work beyond one loop ? 
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We study the field theories for pinned elastic systems at equilibrium and at depinning. Their 
/3-functions differ to two loops by novel "anomalous" terms. At equilibrium we find a roughness £ = 
0.20829804e + 0.006858e 2 (random bond), £ = e/3 (random field). At depinning we prove two-loop 
renormalizability and that random field attracts shorter range disorder. We find £ = |(l+0.14331e), 
e = 4 — d, in violation of the conjecture £ = e/3, solving the discrepancy with simulations. For long 
range elasticity £ = |(1 + 0.39735e), e — 2 — d, much closer to the experimental value (« 0.5 both 
for liquid helium contact line depinning and slow crack fronts) than the standard prediction 1/3. 



The aim of this Letter is to report progress on a concep- 
tual issue and, as a byproduct, to resolve a longstanding 
discrepancy between theory and numerical simulations 
or experiments. The issue is whether it is possible to 
construct a field theory of disordered elastic systems, at 
equilibrium and at depinning, renormalizable beyond one 
loop as for standard critical phenomena. A discrepancy 
exists at present between the value for the roughness ex- 
ponent £ predicted by theory (£ = e/3 exactly) and sim- 
ulations as well as experiments on wetting and on cracks. 

Numerous experimental systems can indeed be mod- 
elled as elastic objects pinned by random impurities, with 
specific features. Interfaces in magnets jjj experience 
either random bond RB (i.e. short range) disorder or 
random field RF (i.e. long range) disorder. Charge den- 
sity waves (CDW) or the Bragg glass in superconduc- 
tors H are periodic objects. The contact line of liq- 
uid helium meniscus on a rough substrate is governed 
by long range elasticity and so are slowly propagating 
cracks §-§. They can all be parametrized by a height 
(or displacement) field u(x) (x being the ci-dimensional 
internal coordinate of the elastic object), with in some 
cases N > 1 components. The roughness exponent £: 



\u{x) - u{x')\ 2 ~ \x-x'\^ 



(1) 



is measured in experiments for systems at equilibrium 
(£ eq ) or driven by a force /. Other exponents describe the 
velocity near the depinning threshold f c , v ~ (/ — f c )^, 
the scaling of the dynamical response, t ~ x z , and the 
local velocity correlation length £ ~ (/ — f c )~ v - 

The study of pinned elastic systems, among a broader 
class of disordered models (e.g. random field spin mod- 
els) , is notably difficult due to dimensional reduction DR 
which renders naive perturbation theory useless ^j,f?]] . In- 
deed to any order in the disorder at zero temperature 
T = 0, any physical observable is found to be identi- 
cal to its (trivial) average in a Gaussian random force 
(Larkin) model. A bold way out of this puzzle was pro- 
posed by Fisher J8| within a one-loop renormalization 
group analysis of the interface problem in d = 4 — e. He 
noted that the coarse grained disorder correlator becomes 



non-analytic beyond the Larkin scale L c , yielding large 
scale results distinct from naive perturbation theory. An 
infinite set of operators become relevant in d < 4, pa- 
rameterized by the second cumulant R(u) of the random 
potential, i.e. V(x, u)V(x', u') — 5 x — x iR(v, — u'). Explicit 
solution of the one-loop Functional RG equation (FRG) 
for R(u) gives several non trivial attractive fixed points 
(FP) to 0(e) proposed in §] to describe RB, RF dis- 
order and in ||, periodic systems (RP) such as CDW 
or vortex lattices. All these FP exhibit a "cusp" singu- 
larity as R*"{u) - i?*"(0)~ \u\ at small Large N 
and variational methods [ll]]^] confirmed the picture and 
the cusp was interpreted in terms of shocks in the renor- 
malized force [Q. A FRG was also developed to one 
loop ||,[l0| to describe the driven dynamics just above 
depinning / = /+, the cusp being linked to the thresh- 
old f c ~ |A'(0 + )|. Surprisingly, the flow equation for 
the correlator A(u) of the force F(x, u) is, to one loop, 
identical to the one of the statics (with A(u) = —R"(u)). 
Extension to temperature T > yielded rounding of the 
cusp in a layer u ~ T and the celebrated creep law p^ ] . 

Despite these successes, serious difficulties remain. 
First, in the last fifteen years since ||, no study has 
addressed whether the FRG yields, beyond one loop, a 
renormalizable field theory able to predict universal re- 
sults HU]. Doubts were even raised |15j about the valid- 
ity of the e-expansion beyond the order 0(e). Second, 
numerous simulations near depinning PJl^ |l8[] seem to 
exclude £ = e/3 argued in |l(| to be exact. In the case 
of long range elasticity, the prediction Q = (2 — d)/3 Q| 
disagrees with the systematically larger value £ « 0.55 
(d = 1) measured for liquid Helium contact line depin- 
ning H and for the in plane roughness of slow crack fronts 
H (see also simulations p9|). 

In this Letter, we address these issues both for dynam- 
ics and statics. The main difficulty is the non-analytic 
nature of the theory (i.e. the fixed point action) at T = 0, 
which makes it a priori quite different from conventional 
critical phenomena. For depinning, we overcome the 
problem and show renormalizability at two-loop order. 
As a result we resolve several questions left unclear in 



1 



previous works. We find that (i) quasi-static driven dy- 
namics differs from statics at two loops (ii) shorter range 
disorder is within the RF universality class and (iii) the 
conjecture £ = e/3 is violated. This last result resolves 
the longstanding discrepancy with simulations. In the 
case of long range elasticity it yields £ « 0.5 for d = 1 
and may thus explain the high value of £ found in ex- 
periments on cracks and wetting. For the statics we find 
ambiguities at T = which we lift using a renormaliz- 
ability condition, yielding fixed points and £ eq to 0(e 2 ). 
This result is also obtained within an independent ex- 
act FRG study j20). The FRG equation for the disor- 
der contains new anomalous terms both for statics and 
dynamics, which are absent in an analytic theory. Our 
predictions for all exponents are shown in Tables 1,11. 
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Table I: exponents for depinning and statics (Ceq) as ob- 
tained, respectively: from setting e = 4 — d in the one loop 
and two loop result, from Pade estimates together with scaling 
relations and from numerical works. For £ eq we have improved 
the estimate using the exact result Ceq(d = 1) — 2/3. 
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Table II: depinning exponents for long range elasticity in 
d = 1: C is consistent with experiments on contact line depin- 
ning (C ~ 0.5 §) and cracks (C » 0.55 ± 0.05 §). 
The starting point is the equation of motion: 



r]d t u x t = dlu xt + F(x, u xt ) 



(2) 



with friction r\ and in the case of long range elasticity 
we replace (in Fourier) q 2 u q by \q\u q in the elastic force. 
Disorder averaged corre lations {A[u x t\) = {A[u x t\)s an d 
responses 5{A[u\) /5h xt — (u xt A[u])s can be computed 
from the standard averaged dynamical action: 



S = 



u xt {r]dt - d x )u xt 



2 

xt "if 

-> u + vt in S. 



Finite temperature is studied adding 
dynamics adding — / J u xt and shifting u 
We study the quasi-static limit v — + , as well as equilib- 
rium dynamics / = where, via fluctuation dissipation 
relations, static quantities can be equivalently computed 



using S or the replicated hamiltonian 21 



It is useful to first study naive perturbation theory, 
in an analytic A(it) i.e. in its derivatives A(™)(0), using 
the diagrammatic rules of Fig. [lj. Since at each vertex 
there are one conservation rule for momentum and two 
for frequency we consider both unsplitted (local a;) and 
splittcd (bilocal t, t') vertices (and splitted a, b vertices in 
the statics). T — power counting yields f t uu ~ x d ~ 2 
and u ~ , where ( = 0(e = 4 — d) has to be deter- 
mined. For an analytic A(it) the perturbation expansion 
of any (analytic) observable yields identical results p^j 
as setting A(u) = A(0) and one obtains the incorrect 
DR roughness £ = e/2. Temperature is formally irrel- 
evant and must be scaled Q as T = f A~ 2+e - 2C with 
the UV cutoff A (and fixed dimensionless T). By power 
counting the only superficially UV divergent irreducible 
vertex functions (IVF) are found to involve only one or 
two response fields u (at T > each T comes with a 
required A 2_d factor to compensate the divergence p3|). 
The statistical tilt symmetry u xt — * u x t + const, (see 
e.g. H[[o]]) further restricts the needed counterterms at 
/ = f c to only one for 77 and one for the full function 
A(u). The one loop (D) and two loops (A,B,C) dia- 
grams which correct the disorder at T = are shown in 
Fig. [l] (unsplitted). The splitted graphs corresponding 
to A in the statics (and which do not vanish or cancel in 
what follows) are shown in Fig. |[ The dynamical dia- 
grams are obtained from the static ones by adding one 
external u on each connected component (e.g. b generates 
61, . . . , be). To escape triviality at T = we must now 
develop perturbation theory in a non-analytic interaction 
A(u) (or R(u)), a non trivial extension of conventional 
field theory. Let us illustrate the new rules. Derivation 
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A(u xt - u xt >) 



FIG. 1. (i) diagrammatic rules for the statics: replica prop- 
agator (« a W(,)o = TSab/q 2 , unsplitted vertex, equivalent split- 
ted vertex — y^ lnh 7^R(u a — ut). (ii) dynamics: response 
= Rq,t-f, unsplitted vertex, splitted vertex 
u xt i) and temperature vertex. Arrows are 
along increasing time. An arbitrary number of lines can enter 
these functional vertices, (iii) unsplitted diagrams to one loop 
D, with inserted counterterm G and two loop A,B,C,E,F. 



propagator {uu)q 
— zUxtUxt' A(u x t 
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FIG. 2. (a-f): the six splitted (static) diagrams correspond- 
ing to two loop A diagram. Below: the corresponding non 
vanishing diagrams in the dynamics. The last one is the only 
non trivial C diagram (see text). 



by extracting a leg from a vertex can be done as usual 
only for a vertex evaluated at a generic u (e.g. graphs 
hi in Fig ||). If it is evaluated at u — (e.g. graph 
ei), one must expand A(u) in powers of \u\, i.e. A(u) = 
A(0) + A'(0+)M + A"(0+)u 2 /2 + ... and carefully apply 
Wick's rules. The result is that the above diagrammatic 
rules (Fig. [|J|) can still be used except that the values 
of the diagrams are different. The graphs of Fig. || cor- 
respond to performing four Wick contractions and some 
end up in evaluating non trivial averages of e.g. sign or 
delta functions. For instance e\, which vanishes in the 
analytic theory since A'(0) = 0, now reads: 

e l = A'(0 + ) 2 A"(u) / R riM Rrid2Rr 3 -riteRr- 3 ,t i F ri ,ti , 
Jti>0,ri 

where F ruU = (sgn(X)sgn(Y)), X = u ri; _ t3 -u ri) _ t4 _ 4l , 
Y = Mo.-t 4 — u o.-t 3 -t 2 i computed with Gaussian aver- 
ages. The limit T -> at v = yields (sgn(X)sgn(F)) = 
fasm((Xy)/ v /(X 2 )(F 2 ))), and a complicated T = 
expression for e\ in the statics p4|]. The opposite 
limit v — > at T — corresponds to depinning, with 
(sgn(X)sgn(F)) -> sgn(t 4 -Ki - £ 3 )sgn(£ 3 +t 2 - U), and 
more generally to A 1 ™) (u t — uy) — > A(™)(v(t — t')) in any 
vertex evaluated at u = 0. 

We now focus on depinning at T = 0. Using these rules 
we compute in perturbation of A = A(u) the contribu- 
tions to the disorder IVF to one and two loops: 

5 1 A = -(A' 2 + (A-A(0))A")I (3) 
S 2 A = ((A - A(0))A' 2 )'7 A (4) 

+i((A - A(0)) 2 A")"/ 2 + A'(0+) 2 A"(/ A - I 2 ) (5) 

with I = J q 1/g 4 and I A = l/gM(?i + q 2 ) 2 [[§, 

whose divergent parts S div A, S div A yield the one loop and 
two loop counterterms respectively. These are computed 
here adding a mass q 2 — > q 2 +m 2 , using dimensional regu- 
larization Im e = N d {\+0(e)), I A m 2e = N d (^ + ±) and 
absorbing N d = (d-2)/(47r) d / 2 r(f ) in A. (|) comes from 
fli + a2 + X)i ^ii t ne fi rs t term in (^J) from all graphs C (not 
detailed) except graph i\ (shown) which contributes to 
the last (anomalous) term in (pf), together with ei, fx, c\ 



(the B contribution vanishes). Inverting the relation be- 
tween bare and renormalized disorders yields the (3 func- 
tion (3 A = dA = eA + eS^A + e{26 2 lv A - 5 l ^A) where 
the 1/e terms cancel nicely, the hallmark of a renormaliz- 
able theory (J 1,1 A is the counterterm to graph G in Fig.l 
and d = —md m ). We obtain the 2-loop FRG equation: 

8A(u) = (e- 20A(u)+(uA>(u)~ [(A( u )-A(0)) 2 ] " 

+~ [(A{u)-A{0))A'{u) 2 }" + l -A'{Q+) 2 A"{u). (6) 

Computing the other needed counterterm, i.e. the renor- 
malized friction r/R = Z~ lr q$, we obtain the dynamical 
exponent z = 2 — dlnZ. The 1/e divergences again can- 
cel yielding the finite result z = 2- A"(0+) + A"(0+) 2 + 
A w (0 + )A'(Q+)(J- ln2). We stress that (|) cannot be 
read at u = p8[. Indeed, it (and the cancellation of 
divergent parts) was established only for u ^ 0. To 
complete two loop renormalizability we checked that IVF 
which are u — quantities are also rendered finite by the 
above counterterms. We found that the time dependence 
in diagrams cancels by subsets as in p2] |, i.e. correlations 
(already rendered finite by the above procedure) are thus 
static for v = + at variance with previous works ^ . 

For periodic A(u) (CDW depinning \U$X§) we find 
a fixed point of (j^) with £ = reading (for a period 
1) A*(u) = 3! + im - (f + ^Hl- u) (0 < u < 1). 
This yields the correlations (u x — uq) 2 = Ad\n\x\ with 
Ad = e/18 + 5e 2 /108, the RP dynamical exponent z = 
2 — -|e — |e 2 and f3 = z/2 from the scaling relation [p|jic|] 

[3 = (z — C)/(2 — C)' Jo 1 ^* becomes non zero to two 
loops, a signature of nonequilibrium effects. 

Another single FP is found to describe both random 
field and all shorter range disorder, including RB, demon- 
strating the instability of the apparent one loop short 
range fixed points. It is determined numerically p4[ ] but 
C is obtained analytically. Integrating (||) over u > 
yields 8D = (e - 3Q D - A'(0+) 3 where D = J + °° A 
(only assuming A(+oo) = 0). The FP condition then 
implies [g (both for RB and RF): 

C = i e + C2£ ^l( 1 + _L ?) = i {1 + 0.14331,), (7) 

where we used that at one loop D* = v / 6e7A*(0) 3 / 2 with 
7 = dy^/y - 1 - In y = 0.54822 This demon- 

strates a violation of the conjecture of [|10[. It reconciles 
theory and numerical results as shown in Table I where 
the dynamical exponent z — 2— |e + e 2 (^ — ^ — j^) = 
2 — |e — 0.04321e 2 as well as (3 obtained via the scaling 
relation, j3 — 1 — — 0.040123e 2 , are also given. 

The case of long range elasticity is obtained chang- 
ing q 2 + m 2 — > \/ q 2 + m 2 in all propagators, shifting the 
upper critical dimension to d uc — 2. It yields a renormal- 
izablc theory, with e = 2 — d and a two loop beta function 
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p4| obtained by multiplying all C(A 3 ) terms in (||) by 
4 In 2. This yields C = 1(1 + 



i.e. 



3^ 1 9 7 \/2 

a strong deviation from e/3 



= 1 



4 In 2 



3 (1 + 0.39735e), 
see Table II), and 



"+ 2 ° 8 ln2 ) = l-fe-0.1133e 2 



9" ' " ^27 7 %/2 108 

We now turn to the statics, using replicas. In the T = 
limit, the FRG beta function at which we arrive 1241: 



dR=(e- 4C oq )i? + Co q uR' 

R"{0))R"' 2 - XR"'{0 



^R" 2 - R"(0)R" 



-\{R" 



! R" 



(8) 



has a new "anomalous" term cx A. The other part, i.e. 
(0) with A = (from graphs a, b and repeated one loop 
counterterm - B graphs cancel in the sum) could as well 
be obtained for an analytic R(u), as in jl4|, which by 
itself would be inconsistent since the FP is non-analytic. 
Ambiguities arise only at two loops (not at one loop since 
R"(0) = R"(0 + )), in the graphs e, / in Fig. | which cor- 
rect R(u) determining A, since some vertices are evalu- 
ated at u = 0. However we have shown that the theory 
can be renormalizable in the usual sense only if: 

A = 1/2 . (9) 

Indeed, the form of the repeated one loop counterterm 
(i.e. to G in Fig.l) S^R = [(R" - R"{0))R"' 2 + (R" ~ 
R"(0)) 2 R"" - R"'{0+) 2 R"]I 2 which is non ambiguous 
because S 1 R(u) is twice differentiable at u = 0, im- 
poses the coefficient of the ambiguous term e + / of 
5 2 R implying (||). Interestingly, this value of A is also 
the only one which prevents the occurrence of a further 
problem in the two loop FRG, the supercusp pS| ]. In- 
deed, e.g. in the periodic case, the FP of (|8|) is R*(u) — 

const. — (j2 + Tos) u2 (l - u ) 2 + J32 _ _ u ) an d 

possesses a stronger singularity than at one loop, since 
R*' is discontinuous. Thus, unless A = 1/2, one has 
J R" = 2R'(0 + ) ^ 0, i.e. a violation of potentiality (as 
naturally occurs above in the driven dynamics). The A = 
1/2 theory yields Ad = jg + j^g for one component Bragg 
glass (and J* A*=0 as natural), £ cq = e/3 for RF disor- 
der and, via numerics, Ccq = 0.20829804e + 0.006858e 2 for 
RB disorder. The corresponding extrapolations (Table I) 
improve the predictions compared to the one loop result. 

Methods aiming at deriving a FRG equation, i.e. com- 
puting A, beyond (physical) renormalizability or poten- 
tiality requirements, are explored in |24f| . An alternative 
exact FRG method, based on multilocal expansion, also 
provides [^0| a procedure to lift the u = vertex ambigu- 
ities at T = 0, and yields (||) with A = _l/2 and universal 
coefficients. A = 1/2 is also recovered |g(J at T > where 
it is easy to see how, at large scale where the running 
temperature Ti flows to 0, anomalous terms as in (^) are 
generated, e.g. from a graph E of Fig. |l| (proportional to 
TiR""(0)R"(u)) since the thermal boundary layer analy- 
sis at one loop @ yields f/i?""(0) -» i?"'(0+) 2 . 

In summary, by finding a way to cope with the diffi- 
culties related to non-analyticity at T = in the FRG, 



we obtained the exponents characterizing depinning and 
statics of pinned elastic systems to next order in e = A—d. 
We predict that similar anomalous terms arise in other 
disordered systems where dimensional reduction fails, 
e.g. random field spin models. 
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At equilibrium, time persistent, i.e. oc 6(ui) quantities 
identify (to all orders in perturbation) to distinct replica 
averages and equal time to same replica averages, 
noting |2(J that graphs can be grouped in subsets (e.g. 
pairs ac, bd, ef in Fig. ^) which vanish by shifting the 
endpoint of an internal line within a splitted vertex, 
similarly the continuum limit in <fi 4 theory requires scal- 
ing as e.g. (j> 6 /A 2 (we thank E. Brezin for this remark). 
P. Chauve, P. Le Doussal and K. Wiese, to be published, 
anomalous contributions (e.g. oc s2) are found to cancel. 
The FP is however unstable in one direction. 
A. A. Middleton Phys. Rev. E 52 R3337 (1995). 
This is why the argument of Ref. |h| (£ = e/3 to all or- 
ders), fails as it misses a contribution in f^°° dudA(u) 
2dD. By contrast at equilibrium, £ oq = e/3 holds to all 
orders if R*'(+oo) 7^ 0, as it should for RF. 
Fractional powers of disorder then appear in equilibrium 
expectation values. 
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